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We formulate the self-consistent separable random-phase-approximation (SRPA) method and 
specify it for Skyrme forces with pairing for the case of axially symmetric deformed nuclei. The 
factorization of the residual interaction allows to avoid diagonalization of high-rank RPA matrices, 
which dramatically reduces the computational expense. This advantage is crucial for the systems 
with a huge configuration space, first of all for deformed nuclei. SRPA takes self-consistently into 
account the contributions of both time-even and time-odd Skyrme terms as well as of the Coulomb 
force and pairing. The method is implemented to description of isovector El and isoscalar E2 gi- 
ant resonances in a representative set of deformed nuclei: 154 Sm, 238 U, and 254 No. Four different 
Skyrme parameterizations (SkT6, SkM*, SLy6, and SkI3) are employed to explore dependence of 
the strength distributions on some basic characteristics of the Skyrme functional and nuclear mat- 
ter. In particular, we discuss the role of isoscalar and isovector effective masses and their relation 
to time-odd contributions. High sensitivity of the right flank of El resonance to different Skyrme 
forces and the related artificial structure effects are analyzed. 



I. INTRODUCTION 

Self-consistent mean-field models with effective energy- 
density functionals (Skyrme-Hartree-Fock, Gogny, rela- 
tivistic) are established as reliable tools for description 
of nuclear structure and dynamics, for a comprehensive 
review see In particular, there is a rising implemen- 
tation of these models to dynamical features of nuclei 
(see, e.g. 0, 0, 0, @), which is caused, to a large ex- 
tend, by exploiting the spectra and reaction rates of ex- 
otic nuclei in astrophysics However, applications of 
self-consistent models to nuclear dynamics are still lim- 
ited even in the linear regime which is usually treated 
within the random-phase-approximation (RPA) . The cal- 
culations are plagued by dealing with high-rank RPA 
matrices. This is especially painful for deformed systems 
where lack of symmetry requires a huge one-particle-one- 
hole {lph) configuration space. 

The RPA problem becomes much simpler if the resid- 
ual two-body interaction is factorized (i.e. reduced to a 
separable form): 

K 

(pxhi\V res \ph) a+ a^apaft — > ^ i^k,k'X k X k > 

where X k = J2 P h(p\Xk\h)apO-h is a hermitian lph op- 
erator and Kkk' is a matrix of strength constants. The 
factorization allows to reduce a high-rank RPA matrix 
matrix to much smaller one with a rank 4K (where the 
coefficient 4 is caused by the isospin and time-parity fac- 
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tors, see discussion at the next Section). The separa- 
ble expansion can be formulated in a fully self-consistent 
manner and provide a high accuracy with a small number 
(K ~ 2 — 6) of the separable terms. 

Several self-consistent schemes for separable expan- 
sions have been proposed during the last decades H 0, 
OH HE UA EE LL-1 Ull • However, these schemes were not 
sufficiently general. Some of them are limited to ana- 
lytical or simple numerical estimat es [1 , 0, ITol . [TT| , the 
others are not fully self-consistent |13j . Between them 
the approach |14L for Skyrme forces is quite promis- 
ing. However, it still deals with RPA matrices of rather 
high rank (~ 400). 

Recently, we have developed a general self-consistent 
separable RPA (SRPA) approach applicable to arbitr ary 
density- and current-dependent functionals 0, 0, 0, 
[H |M El E2. The method was implemented to the 
Skyrme functional Q, I2H . I24I for both spherical 
OH [H HJ EH and deformed [12 nuclei. In SRPA, 
the one-body operators X k and associated strengths K k , k > 
are unambiguously derived from the given energy-density 
functional. There remain no further adjustable parame- 
ters. The success of the expansion depends on an appro- 
priate choice of the basic operators X k . Experience with 
spherical SRPA gives guidelines for an efficient choice 
such that a few separable terms suffice to reproduce accu- 
rately the exact RPA spectra [la 13123 ■ The success be- 
comes possible due to the following factors: i) an efficient 
self-consistent procedure 0, 0] based on sound physical 
arguments, ii) proper inclusion of all parts of the resid- 
ual interaction, time-even as well as time-odd couplings, 
iii) incorporation of the symmetries (translation, particle 
number, ...) leading to a correct description of the related 
zero-energy modes, iv) building the separable operators 
by such a way that they have maxima at different slices 
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of the nucleus and thus cover both surface and interior 
dynamics. Furthermore, we note that SRPA is very gen- 
eral and can be applied to a wide variety of finite Fermion 
systems. For example, SRPA was derived for the Kohn- 
Sham functional |2y, |2jj and widely used for description 
of linear dynamics of valence electrons in s phe rical and 
deformed atomic clusters E [H IH H US Ml ■ 

The enormous reduction of computational expense by 
SRPA is particularly advantageous for deformed systems 
where lph configuration space grows huge. SRPA be- 
comes here a promising tool for large scale studies. It 
is the aim of this paper to present a first implementa- 
tion of SRPA for axially deformed nuclei with pairing. 
Like the full RPA, SRPA follows two strategies to com- 
pute the dynamical response. It can be calculated via 
determination of RPA eigenvalues and eigenstates or by 
a direct computation of the multipole strength functions 
related to experimental cross sections. We will discuss 
both ways. 

As a first test, we will apply the method to descrip- 
tion of isovector El and isoscalar E2 giant resonance 
(GR) in the deformed nuclei 154 Sm, 238 U, and 254 No. 
These nuclei cover a broad size range from rare-earth 
to super-heavy elements. Four different Skyrme forces 
(SkT6 [13, SkM* [H, SLy6 0, and SkI3 |a3) will be 



used to explore dependence of the results on the actual 
parametrization. For our aims it is important that these 
forces represent different values of some relevant nuclear 
matter characteristics (isoscalar and isovector effective 
masses and asymmetry energy). We will discuss depen- 
dence of the collective strength on these characteristics, 
scrutinize the impact of time-odd coupling terms in the 
Skyrme residual interaction and demonstrate important 
role of the Landau fragmentation. 

The paper is organized as follows. The general SRPA 
formalism is presented in Sec. Inland specified for Skyrme 
forces in Sec. IIIII In Sec. HVI we discuss the choice of the 
input operators. Results of the calculations are analyzed 
in Sec. [3 The summary is done in Sec. I VII Some 
important details of the method are given in Appendices 
A-C. For more details, see documentation in |22| . 



II. BASIC SRPA EQUATIONS 

A. The separable expansion 

In the most general form, the factorization of the resid- 
ual interaction (Q) reads 



ss' k,k' — l 



'X s kX s 'k' + TJsktS'k'YskYs'k'} 



(2) 



where the indices s and s' label neutrons and protons, 
X s k are time-even hermitian one-body operators and Y sk 
are their time-odd counterparts, n s k,s'k' and rj^ ^^i are 
the strength matrices. Time reversal properties of the 
operators are 

TX sk T- 1 = j*X sk , 7 # = +l , (3) 
TY sk T- 1 = rfy sk , -yjf = -l , 

where T is the operator of time reversal. The expan- 
sion © needs to take care of both classes of operators 
since the relevant Skyrme functionals involve both time- 
even and time-odd couplings, see Q and Appendix f£] 
Though time-odd variables do not contribute to the static 
mean field Hamiltonian of spin-saturated systems, they 
can play a role in time-dependent perturbations and nu- 
clear dynamics. Altogether, the presence of time-even 
and time-odd couplings in the functional naturally leads 
to a formulation of the separable model in terms of hermi- 
tian operators with a given time-parity. These operators 
have the useful property 

([!,£]} ~(l- 7 #7#) (4) 



which means that the average value of the commutator 
at the ground state |) is not zero only for operators of 
the opposite T-parities (7^ = —7^). This property will 
be widely used in the following. 

B. Linearized time-dependent mean field 

RPA is the limit of small-amplitude harmonic vibra- 
tions around the ground state. The dynamics is formu- 
lated in general by a time-dependent variation on the 
basis of a given energy functional E(Jf(r,t)) 

(*(t)|£|tf(t)) — * E(J?(f,t)) = J H(r,t)dr. (5) 

We deal with a set of densities Jf (r) where (3 denotes the 
type of density (spatial density, kinetic density, current, 
spin density, spin-orbit density, ...) and s labels protons 
and neutrons. For reasons of compact notation, we com- 
bine the density type (3 with the space point r into one 
index a such that 

a =09,0 , £... = W<F... (6) 

a p J 
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The densities are related to the corresponding one-body 
operators J" (see the list in Appendix [X| as 



Jf(t) = (*(t)|J*|*(t)> 



(7) 



Further, the state ^(t)) is the underlying Bardeen- 
Cooper-Schrieffer (BCS) mean-field state composed from 
the single-particle states ip h (r , t) and the corresponding 
pairing occupation amplitudes v^. The time-evolution is 
determined by variation of the action (t)\idt\^> it)) — 
J drH(f, t). Till now we keep the occupation amplitudes 
Vh fixed at their ground state values and consider the 
variation of the single particle states. This yields the 
time-dependent mean-field equations as 



at 



(8) 



with the mean field Hamiltonain h being a functional of 
the local and instantaneous densities Jf(r, t). The freez- 
ing of the occupation amplitudes vh somewhat inhibits 
application of SRPA for vibrational modes with a strong 
pairing impact (e.g. for low-lying modes in neutron-rich 
light deformed nuclei) |3(j- However, in the present study 
we are interested in giant resonances where pairing dy- 
namics plays a minor role. 

In the linear regime of small amplitude oscillations, the 
time-dependent state consists of the static ground state 
|) and a small time-dependent perturbation 



|*(t)> = |) + |**(t)> 



(9) 



where both |\l/(t)) and |) are BCS states. Hence, all dy- 
namical quantities can be decomposed as a sum of the 
stationary ground state and small time-dependent parts: 



related to the oscillations of the system. Note that the 
second functional derivative is to be taken at the ground 
state density as indicated by the index J = J. This 
holds for all second functional derivatives and so we will 
skip this explicit index in the following. For the brevity 
of notation, we will also skip the explicit dependencies 
on space coordinates and come back to these details in 
Sec. IIIII where the residual interaction for the Skyrme 
functional is worked out. 

The linearized equation of motion reads 



(15) 



Further steps require a more specific view of \5^!). This 
will be done in the next subsections from different points 
of view, macroscopic and microscopic. 



C. Macroscopic part of SRPA 

1. Scaling perturbed wave function 

It is convenient to obtain the perturbed mean-field 
state |W(f)) by the scaling transformation ^(| 



!*(*)>. 



K 

n exp t 

fe=i 



-iq s k(t)P s k] exp[-ip s k(t)Q s k]\)s (16) 



where both | s and \) s are the Slater determinants 
and Q s k (r) and P s k (f) are generalized coordinate (time- 
even) and momentum (time-odd) hermitian one-body op- 
erators. These operators fulfill the properties 



Q sk = Qt 
P sk = i[H,Qsk}= P s 



sh- • 



1 Q T 



1. 



(17a) 
(17b) 



J°(r,t) = J?(f)+5J?(f,t), (10a) 
SJ«(r,t) = (*(t)|#|*(t)> - (J?). (10b) 

Inserting l|10[l into (J5J and keeping the terms up to the 
first order in 8Jf(r,t), one obtains the single-particle 
Hamiltonian 



h{t) = h + h rcs (t) 
with the static mean-field part 



h 



Y—J a 



and the time-dependent response 



E 



5 2 E 



5J?5J5' 



(11) 



(12) 



j=j 



(13) 

J?6J2'(t) (14) 



where H = ho + V ICS stands for the full Hamiltonian em- 
bracing both the one-body mean-field Hamiltonian and 
the two-body residual interaction. The commutator in 
(|17b(l is assumed to be mapped into the one-body do- 
main (see, e.g., the mapping into h ros in Eq. (11911 1. If 
the functional includes only time-even densities, then V res 
does not contribute to the commutator and so H can be 
safely replaced by ho. 



2. Separable operators and strength constants 

The operators l|17fl generate time-even and time-odd 
real collective deformations q s k(t) and p s k(t)- Using Hlbj) 
and assuming only small deformations, the transition 
densities (110 bl) read 



k 

+ P sk(t)([QskJ:})} 



(18) 
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where, following (J2J, time-even densities contribute only 
to responses (|[P s fc, ]|) while time-odd ones only to 
(|[Qsfc> ^"]|)- Then the response Hamiltonian lj*H5|) re- 
casts as 



he S (t) = ^2{q S k(t)x sk +p sk (t)Y sk } 



(19) 



sk 



where the time dependence is concentrated in the ampli- 
tudes q sk {t) and p sk (t) while all time-independent terms 
are collected in the hermitian one-body operators 



X s k 



i V JJL-([P ak ,j?])J«', (20a) 



<H, SJZ'SJ? 



with the properties 

X = x\ 
Y = Y\ 



7? = 



+1, 
-1, 



X* = X, 
Y* = -Y. 



(21a) 
(21b) 



SJ5SJ: 



Obviously X sk and Y sk are the reasonable candidates 
for the time-even and time-odd operators in the separable 
expansion J2J. The operator X sk involves contributions 
only from the time-even densities while the operator Y sk 
only from the time-odd ones. The upper index s' in the 
operators 11201) determines the isospin (proton or neutron) 
subspace where these operators act. This is the domain 
of the involved density operators J" . 

To complete the construction of the separable expan- 
sion (J2J, we should yet determine the matrices of the 
strength constants n s k,s'k' and Tj s k, s 'k'- This can be done 
through variations of the basic operators 



5X sk {t) = ^>{t)\X sk \^{t)) - (X sk ) = iJ2<ls>k'(t)([P s , k ,,X s s ' k ]) = -$> s ' fe '(*K~it', 



sk s 



s'k' 



s'k' 



8Y sk {t) = (*(t)|y >fc |*(t)) - (Y sk ) =iJ2Ps'k'(t)([Q s 'k',Y s i}) =-Dm'(*)4 



sk 



s'k' 



s'k' 



(22a) 
(22b) 



where 



s'k'sk sk,s'k' 



^s'k'.sk ^sk 7 s'k' 



i{[p s , k ,,xi]) = Y,jj^^ < 

aa' s s 
aa' s s 



(23a) 
(23b) 



Eqs. 1(23(1 represent the elements of the symmetric ma- 
trices which are inverse to the matrices of the strength 
constants in (j2J. Indeed, Eqs. (|22(l can be recast to 

-}^K s 'k',skSX sk (t) = q s 'k'(t) , (24a) 

S fe 

-£f7 a 'k',«*<* y «fc(*) = Ps'k'(t) ■ (24b) 

From that we read off the response Hamiltonian 1191) as 

h rcs (t) = - ^ ^ {K s 'k',sk$X sk (t)X s , k < (25) 

s'k' sk 

+Vs'k' ,sk^Y sk (t)Y s r k 'j , 

which leads within the collective space of the generators 
{Psk i Qsk} to the same eigenvalue problem as the sep- 
arable Hamiltonian 

H = h + V r fP (26) 
with V^f given in Eq. @ (see dEl). 



In principle, we already have in our disposal the macro- 
scopic SRPA formalism for linear regime of the collective 
motion in terms of real harmonic variables 

q.h(t) = q sk cos(uit) = ^q sk {e lut + e~ wt ) , (27a) 

Psk(t) = Psk sin(wt) = ±:Psk(e iwt - e" Mt ) . (27b) 

Indeed, Eqs. I|20|l and 1)23(1 deliver the one-body opera- 
tors and strength matrices for the separable expansion of 
the two-body interaction. By substituting the response 
Hamiltonian (|25|l and the perturbed wave function l|l(j|) 
into time-dependent HF equation 1(15(1 one gets the eigen- 
value problem. The number K of the collective variables 
(and thus of the separable terms) is dictated by the accu- 
racy we need in the description of collective modes. For 
K = 1, the method in fact is reduced to the sum rule 
approach with one collective mode (38|. For K > 1, we 
have a system of K coupled oscillators and the method 
becomes similar to so-called local RPA j3^, suitable 
for description of branching and gross-structure of collec- 
tive modes. 
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D. Microscopic part of SRPA 

The macroscopic SRPA as outlined in section III CI 
serves here as a convenient tool to derive the optimal 
separable expansion. But it cannot describe the Landau 
fragmentation of the collective strength. For this aim, we 
should build the microscopic part of the model. In what 
follows, we will consider the eigenvalue problem and the 
direct computation of the strength function. 

The perturbation \S^f) belongs to the tangential space 
of the variations of a mean field state. For the pure 
Slater states, they are all conceivable one-particle-one- 
hole (lph) excitations 



A p h — a l a h 



A P h — ajjflp 



(28) 



In the BCS case, the elementary modes are reduced to 
the two-quasiparticle (2qp) excitations A ph 



A] 



;,Ut 



a-aj where a] generate the BCS quasiparticle state j 

and j is its time reversal. Just this case is employed 
in our actual calculations. Both lph and 2qp excitations 
have much in common and result in the same microscopic 
SRPA equations (for exception of pairing peculiarities 
outlined in Appendix Q . So, in what follows, we will 
not distinguish these two cases. In particular, we will 
use for both excitations one and the same notation ph. 



1. Eigenvalue problem 

To formulate the eigenvalue problem, we will exploit 
the standard RPA technique for the separable Hamil- 
tonian i|26|) where the separable operators and strength 
constants are delivered by the macroscopic SRPA, see 
Eqs. and J23). Following 0, the collective motion 
is represented in terms of 



!*(*)>* 

cl 



exp 



C v e 



-\-iu) u t 



|), (29a) 



EE te^-W ' (29b) 



s ph£s 



where Cl creates the one-phonon eigenmode v\ the oper- 
ators A' h and A ph are defined in i|28|): the perturbed 
|^(t))i, and ground |) states have the form of Slater 
determinants. We employ here the Thouless theorem 
[40| which establishes connection between two arbitrary 
Slater determinants. The wave function (|29|l is a mi- 
croscopic counterpart of macroscopic scaling ansatz (1 1 01) . 
But here we aim a fully microscopic description of exci- 
tations covering all the lph space while in the previous 
sections we considered macroscopic flow as a benchmark 
for forming the separable interaction. 



The time-dependent response Hamiltonian h res (t) for 
the mode C\,e~ %U3vt and the separable interaction (J2J 
reads 



<r sk 



+E^ fc E^fc'([n'fc',^]> , (so) 



sk s' k' 



where 



Isk 



^ ] K sk.s'k' 



(31a) 



s' k' 



£ (ph\x!: k ,)( c ;-+ c ;+), 

s" phGs" 

Psk = i^2vsk,s<k< (31b) 

s'k' 

' " v— v-\ 



E E (p*i*&>(« 



ph ^ph 



s" ph£s" 



Note that here the values q" k and p v sk are new objects 
generated by the eigenmode C^. They serve for notation 
of /i res and have no relation to the collective generators 
qsk, Psk used in the section III Cl Substituting the ansatz 
B29JI into the linearized time-dependent HF equation l|15|) 
and employing the form H3()|) for the response field yields 
the expression for lph expansion coefficients 



^phes 



E 

s'k' 



2(e p h ± w v ) 



, (32) 



where e p h is the unperturbed energy of lph pair. In the 
derivation above we used the operator properties 121|l . 
Then the matrix elements (ph\X^ k ,) and {ph\Y5,,) are 
real and image, respectively, while all the unknowns c"^, 

Cfc and Psk are re- 
inserting the result (|32|l into Eqs. (|31(l yields finally 
the set of SRPA equations for the values q" k an d P v sk 

E\q v _\p( xx \ _ K ~} 1 _i_ f) v -F^ XY ^- \ = 
l^sfcL s' k' ,sk sk,s'k'* "sk s'k',sk> ' 

sk 

E ^sk F i'k',lk + P V ski F{ s'k',lk ~ ^skyk'}} = (33) 

sk 

where 
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F ^ k B ,lk = "AB E E 



1 



a phGs P h 



{(ph\Ai- k )*(ph\B s a , k ,)(e ph + w v ) + (ph\M- k )(B s s>k ,\ P h)(e ph - uj v )} (34) 



r 



with A, B G X, Y and 



CtAB 



for A 
for A ■ 
for A - 




The system of linear homogeneous equations l|33|l has a 
non-trivial solution only if its determinant is zero. This 
yields the dispersion equation to obtain the RPA eigen- 
values LU,,. 



2. Strength function 



states |4l|. For this aim, we will use the Cauchy residue 
theorem. Namely, the strength function is represented as 
a sum of the residues for the poles z — ±w„ . Since the 
sum of all the residues (covering all the poles) is zero, 
the residues with z — ±io v (whose calculation is time 
consuming) can be replaced by the sum of residues with 
z — lj ± z(A/2) and z = ±s p h whose calculation is much 
less expensive. The explicit derivation is given in |22| . 
The final expression reads 



z L J2f30>F w/ (z)A p (z)A p ,(z) 



nF(z) 



- z=u>HA/2 



When exploring the response of the system to time- 
dependent external fields, we are often interested in the 
total strength function rather than in the particular RPA 
states. For example, giant resonances in heavy nuclei are 
formed by thousands of RPA states whose contributions 
in any case cannot be resolved experimentally. In this 
case, it is more efficient to consider a direct computation 
of the strength function which avoids the details and cru- 
cially simplifies the calculations. 

For an external electric field of multipolarity EX/i, we 
define the strength function as 

S l (X^lu) =Y J "iMl llv t(u ~ (35) 



where 



C(cj - U) v ) 



A 



(36) 



2tt {uo - lu,) 2 + (A/2) 2 
is the Lorentz weight with an averaging parameter A and 



M x „ v = — 



1 

V2 



ph£s 



is the matrix element of E\\i transition from the ground 
state to the RPA state \v). Further, e e J^ is an effective 
charge to be specified later. The operator of the electric 
external field in the long-wave approximation reads 



1 



(38) 



The (normalized) amplitudes follow from l|32|l . Un 



like the standard definition of the strength function with 
5(u) — w„), we exploit here the Lorentz weight. It is con- 
venient to simulate various smoothing effects. 

The strength function l|35|l can be recast to the form 
which does not need information on the particular RPA 



ph£s 



e ph) 



(39) 



where 3 means the image part of the value inside the 
brackets, F(z) is the determinant of the RPA matrix l)33|) 
with uj v replaced by the complex argument z, Fppi(z) is 
the algebraic supplement of the determinant, and 



4?M 



E^ 



E 

ph^s' 



e ph (ph\X s J k ){ph\f Xp ) 

4h - z2 



i J2 e e J f 



E 



Lo v (ph\Y s j){ph\f Xil ) 

4h - z2 



(40a) 



(40b) 



For the sake of brevity, we introduced in l|39|l the new 
index [3 — {skg} where g—1 for time-even and 2 for 



For example, A s k g =i = A 



- aW 



sk 



and 



time-odd quantities. 

The first term in l|39() collects the contributions of the 
residual interaction. It vanishes at V^ cs = 0. The sec- 
ond term is the unperturbed (purely two-quasiparticle) 
strength function. 



E. Basic features of SRPA 

Before proceeding to further specification of the 
method, it is worth to comment some its essential points. 
• To determine the unknowns c"^ s of a full (non- 
separable) RPA, one requires diagonalization of the ma- 
trices of a high rank equal to the size of the lph basis. 
The separable approximation allows to reformulate the 
RPA problem in terms of a few unknowns q$, and p k and 
thus to reduce dramatically the computational effort. As 
is seen from (|33|l , the rank of the SRPA matrix is equal 
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to AK where K is the number of the separable operators. 

• The number of the SRPA states \v > is equal to the 
number of the relevant lph configurations used in the 
calculations. In heavy nuclei, this number can reach 10 4 - 
10 5 . Every state \v > is characterized by the particular 
set of the values q" sk and p v sk . 

• Eqs. J20|) and relate the basic SRPA values with 
the initial functional and input operators Q s k- After 
choosing Q s k, all other SRPA equations are straightfor- 
wardly determined following the steps 



Qsk => ([Qsk, J?]) =>■ Y skl r) s l s , k , 



Xsk, K sk , k ,. 



(41) 



As is discussed in Sec. II VI the proper choice of Q sk is 
crucial to achieve a good convergence of the separable 
expansion with a minimal number of separable oper- 
ators. SRPA itself does not provide a recipe to get Q s k 
but these operators can be introduced following intuitive 
physical arguments. 

• SRPA restores the conservation laws (e.g. translational 
invariance) violated by the static mean field. If a symme- 
try mode has a generator P sym , then we keep the conser- 
vation law [H, -P S ym] = by including P sym into the set 
of the input generators P s k together with its complement 

Qsym — ^[-^\ -fsym] ■ 

• The basic SRPA operators can be expressed via the 
separable residual interaction (J2J as 

X sk = Psk] P h, tk = -*[Kcs P : Qek]ph (42) 

where the index ph means the lph part of the commu- 
tator. It is seen that the time-odd operator P s k retains 
the time-even part of V t s °p to build X s k ■ Vice versa, the 
commutator with the time-even operator Q s k keeps the 
time-odd part of V^F to build Y sk . 



III. SRPA WITH THE SKYRME FUNCTIONAL 
A. Skyrme functional 

We use the Skyrme functional [2j| in the form 0, HM 



E = / dr (7Y kin + 7Y S k (fi s , r a , a s , j s , 3 S ) (43) 
+Hc(Pp) + H pail (x s )) , 



where 



^kin — - — t , 
2m 

e 2 f 1 
Hc = ~2 J d ^Pp^)y—r-\P P {?i) 

- Je 2 (^Mr)]t, 

4 7T 



(44) 

(45) 
(46) 



Hs k = |p 2 -|E^ ( 47 ) 

S 

-^(Ap) + |]>>(Ap s ) 



, ^3 a +2 ^3 a 2 
s 

+ b 1 ( P T-?)~b> 1 Y,(PsTs-j 2 s ) 

s 

-64 (p(VQ) + a ■ (V x J)) 

-K Yl (p>&$.) + a s • (V x j s ) 



arc kinetic, Coulomb, pairing and Skyrme terms respec- 
tively. The densities and currents used in this functional 
are defined in the Appendices El an d O Densities with- 
out the index s involve both neutrons and protons, e.g. 
p = p p + p n - Parameters 6$ and a are fitted to describe 
ground state properties of atomic nuclei, see e.g. 0. 

For the sake of brevity, we omit here derivation of the 
Skyrme mean field which can be found elsewhere, e.g. in 
[TsL I22I l38l | , but present only the main values entering 
the SRPA equations. 



B. Second functional derivatives 

The crucial ingredients of the Skyrme SRPA residual 
interaction are the second functional derivatives enter- 
ing expressions for the basic operators (|20|l and strength 
matrices l|23|) . They read 



S 2 E 



5p Sl {ri)6p s {r) 

(a + 2) (q + l) Q 
+03 ; 



bo - b' S SSl - (6 2 - b' 2 5 SSl )Ap 1 
P a (r) 



, r a(a- 1) a -2/^sr 2 
-h [ o P ( r ) P^ V) 



2a 



p a -\r){p s (r)+p Sl (r)) 



+5 SS1 -P a (7)} 

-SssApli-y^Ppir)}- 213 } 5(r- fi) 



5 2 E 
5T ai (fi)dp a (f f ) 

5 2 E 
5^s Sl (ri)6p s (r) 

6 2 E 
5/o Bl (r=i)(5S s (r) 



|f- fi| 



(48) 



[h - b[d ssl ]S(r ~ n) , (49) 
= [6 4 + & 4 ^ 1 ]V ?1 <5(f-ri) , (50) 

: -[& 4 + fo 4 <W]V r - 1 £(f-r 1 ) (51) 



8 



for time-even densities and 
S 2 E 



2[-bi + b' 1 5 SSl ]5 mn 5(fi ~ r), 
S 2 E 



-[b 4 + b' i S ssl ]e m „i'Vf l! i8(rx - r) 



(52) 



(53) 



for time-odd densities. The last two terms in (|48|) rep- 
resent the exchange and direct Coulomb contributions. 
Further, the indices m, n in H52[) and (|53|l run over the 
three basis spatial directions of the chosen representation 
(in our case the cylindrical coordinate system) and e mn i 
is the totally antisymmetric tensor associated with the 
vector product. 



C. Presentation via matrix elements 

The responses in l|2Uf) and (|23[) are expressed in terms 
of the averaged commutators 



([A,B]} with 7# = -7# 



(54) 



where |) is the quasiparticle vacuum. Calculation of these 
values can be considerably simplified if to express them 
through the matrix elements of the operators A and B. 

To be specific^ we associate the operator A with a time- 
even operator Q or a time-odd operator P, which have 
real or imaginary matrix elements, respectively. Then 

([Qsk,B s ]) = Ai (ph\Q sk )%(ph\B s ) , (55a) 

ph<E& 

{[P sk ,B s ]) = -4 (ph\P sk )$i(ph\B s ) , (55b) 

ph£s 

where both p and h run over all single-particle states 
and 9 and 3? mean the imaginary and real parts of the 
values to the right. The pairing factors are included in 
the single-particle matrix elements. 

Then elements of the inverse strength matrices are real 
and read 



ph£s' 



(56a) 



%<\>,sk = -i([Q*,Y k ) = 4 WQs'k') %(ph\Y s i}) . 

ph£s' 

(56b) 

The responses entering X and Y in (12011 are also real and 
read 

K a x , sk = i([p sk , = -a & h \ p °») a<pfci#> , 

ph^zs 

(57a) 



n^ sk = i{[Q sk , = -*Yl (p*I&*> ^(p h \ J s) 

ph£s 

(57b) 

where (ph\Jf\) are transition densities. 

We actually deal with axially symmetric systems. The 
modes then can be sorted into channels with a given 
angular-momentum projection p to the symmetry axis z. 
For even-even nuclei p takes integer values. The explicit 
expressions for the responses in cylindrical coordinates 
(see definition of these coordinates in the Appendix lB|) 
then read 

IyMV = i{[QskJ s }) (58) 
= (e P jy :Sfe (p, z) + ezjv,sk(P: z )) cosfiO 

sy,3k(r) = i([Qsk,s s }) (59) 
= (e P s Y ^ k (p,z) + e z s^ sk (p,z))sm f ie 

+ e e s^r sk (p,z) cos fj,8 , 

Px, sk {r) = i([P sk ,p s }) = Px, sk {p, z)cosp9 , (60) 

Tx, s k(r) = i([P s k,T s ]) = T X , sk (p, z)cosp9 , (61) 

3x, sfc (r) = *<[A*A]> (62) 

= (e P ^x,s k (P> z ) + e^x.skiPi z )) cos P 

+ ee%x,s k {P, z ) sin P° 

where {p, z}-depending response components are real. 
All the dependence on the spatial coordinates follows 
from the transition densities entering the responses. 

The response components involving sin p9 obviously 
vanish for modes with p = 0. 

Explicit expressions in cylindrical coordinates for 
strength matrices, responses, transition densities, matrix 
elements, Coulomb contributions and other SRPA values 
can be found in j^. It is to be noted that the present 
calculations do not yet include the Coulomb contribu- 
tion to the residual interaction. This introduces the un- 
certainty reaching up to ~0.4 MeV in an average peak 
position • Such effect is safely below the precision of 
the present investigation. 

IV. CHOICE OF INITIAL OPERATORS 

As was mentioned in Sec. Ill El SRPA starts with the 
choice of appropriate generating operators Q sk , see the 
sequence of the model steps in (|41H . The SRPA formalism 
itself does not provide these operators. At the same time, 
their proper choice is crucial to get good convergence 
of the separable expansion (J2J with a minimal number 
of separable terms. The choice should be simple and 
universal in the sense that it can be applied equally well 
to different modes and excitation channels. 

We propose a choice inspired by physical arguments. 
The main idea is that the generating operators should 
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explore different spatial regions of the nucleus, the sur- 
face as well as the interior. The leading scaling genera- 
tor should have the form of the applied external field in 
the long-wave approximation, which is most sensitive to 
the surface of the system. Since nuclear collective mo- 
tion dominates in the surface region, already this gener- 
ator should provide a good description. Next generators 
should be localized more in the interior to describe an 
interplay of surface and volume vibrations. For EX giant 
resonances in spherical nuclei, we used a set of generators 
with the radial dependencies in the form of power and 
Bessel functions, [lg- In the present study for deformed 
nuclei, we will implement, for the sake of simplicity, only 
generators with the power radial dependence 

4(f) = r A+3 ( fc - 1 ) 0^.(0) + h.c). (63) 

The separable operators and Yj, with k = 1 are mainly 
localized at the nuclear surface while the operators with 
k > 1 allow to touch, at least partly, the nuclear inte- 
rior. This simple set seems to be a good compromise for 
the first calculations of the giant resonances. Our anal- 
ysis shows that already two first operators with k = 1,2 
suffice for a spectral resolution of 2 MeV, as discussed 
below. In the next studies we plan to enlarge the list 
of the input generators so as to cover properly the nu- 
clear interior and to take into account the coupling of the 
modes with different A, induced by the deformation. 



V. RESULTS AND DISCUSSION 

A. Aims and details of calculations 

The main aim of the present calculations is to demon- 
strate ability of SRPA to describe multipole giant reso- 
nances (GR) in deformed nuclei and to explore depen- 
dence of the resonance strength distributions on the nu- 
clear matter properties and some terms of the Skyrme 
functional. The isovector giant dipole resonance (GDR) 
and isoscalar giant quadrupole resonance (GQR) in the 
axially deformed nuclei 154 Sm, 238 U, and 254 No will 
be considered. We will employ a representative set of 
Skyrme forces, SkT6 H2, SkM* H|, SLy6 0, and SkI3 
|35| . with different features (isoscalar and isovector ef- 
fective masses, asymmetry energy, ...). These forces are 
widely used for description of ground state properties and 
dynamics of atomic nuclei [lj including deformed ones. 
As is seen from Table 1, the forces span a variety of nu- 
clear matter properties while all correctly reproducing 
experimental values of the quadrupole moments in 154 Sm 
and 238 U. 

The calculations use a cylindrical coordinate-space grid 
with the spacing of 0.7 fm. Pairing is treated at the BCS 
level and is frozen in the dynamical calculations. The col- 
lective response for the GDR (A = 1) and GQR (A = 2) 
is computed with two input operators (|63[1 at k=l and 2. 
Both GR are calculated in terms of the energy-weighted 



TABLE I: Nuclear matter and deformation properties for the 
Skyrme forces under consideration. The table represents the 
isoscalar effective mass mo/m, symmetry energy a sym , sum 
rule enhancement factor k, isovector effective mass m\/m = 
1/(1 + ft), and quadrupole moments Q2 in 154 Sm, 238 U, and 
254 No. Experimental values of the quadrupole moment in 
154 Sm and 238 U are 6.6 b and 11.1 b, respectively HJ. 



Forces 


mj/m 


a sym [MeV] 




m*/m 


< 


22 [b] 














154 Sm 


238 y 


254 No 


SkT6 


1.00 


30.0 


0.001 


1.00 


6.8 


11.1 


13.7 


SkM* 


0.79 


30.0 


0.531 


0.65 


6.8 


11.1 


14.0 


SLy6 


0.69 


32.0 


0.250 


0.80 


6.8 


11.0 


13.7 


SkI3 


0.58 


34.8 


0.246 


0.80 


6.8 


11.0 


13.7 



(L=l) strength function (|39[) with the averaging parame- 
ter A= 2 MeV (as most suitable for the comparison with 
the experiment). The factorization of the residual inter- 
action and the strength function technique dramatically 
reduce the computational effort. For example, by using 
a PC with CPU Pentium 4 (3.0 GHz) we need about 
25 minutes for the complete calculations of the GDR in 
238 U, including 1 minute for computation of the strength 
function itself. 

The isovector GDR and isoscalar GQR are calculated 
with the effective charges e?JJ = N/A, e%* = —Z/A and 
eyf = ef/f = 1, respectively. The isoscalar dipole spuri- 
ous mode is located at 2-3 MeV. The deviation from the 
desirable zero energy is caused by several reasons. First, 
we have neglected in the present study the contribu- 
tion from the Coulomb residual interaction. The second 
(more important) reason is that the nucleus is treated in 
a finite coordinate box, which artificially binds the cen- 
ter of mass. Larger numerical boxes could help here but 
at quickly increasing computational cost. In any case, it 
suffices for our present purposes that the center-of-mass 
mode lies at a low energy and thus is safely separated 
from the GDR. 

We use a large configuration space including the single- 
particle spectrum from the bottom of the potential well 
up to - +16 MeV. This results in 7000-10000 dipole and 
11000-17000 quadrupole two-quasiparticle configurations 
in the energy interval - 100 MeV. The relevant energy- 
weighted sum rules are exhausted by 85 — 95%. Such a 
basis is certainly enough for the present aims. 



B. Discussion of results 

Results of the calculations are presented in Figs. 1-3. 
The first two figures compare the calculated GDR and 
GQR with the available photoabsorption 0> @ an d 
(a, a') [5(J experimental data. It is seen that all four 
Skyrme forces provide in general an appropriate agree- 
ment with the experiment. So, SRPA indeed can give a 
robust treatment of the multipole GR and, what is im- 
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FIG. 1: The isovector GDR in lj4 Sm, 238 U, and 254 No cal- 
culated with the Skyrme forces SkT6, SkM*, SLy6 and SkI3. 
The plots depict: the collective strength calculated with two 
(solid curve) and one (dotted curve for 254 No) input opera- 
tors, the unperturbed quasiparticle strength (dash curve) and 
the photoabsorption experimental data |48l.l49f (triangles). 

portant, does this with a minimal computational effort. 

To illustrate accuracy of the method, we give for 254 No 
the strength functions calculated with one (k=l) and two 
(k=l,2) input operators. It is seen that both cases are 
about equal for the GQR with its simple one-bump struc- 
ture. At the same time, the second operator considerably 
changes gross structure of the more complicated GDR. 
We have checked that inclusion of more operators (k > 
3) does not result in further significant modification of 
the GDR. So, the approximation of two input operators 
seems to be reasonable, at least for the present study at 
a resolution of 2 MeV. 

The figures also show the unperturbed strengths, i.e. 
the mere two-quasiparticle (2qp) distributions without 
the residual interaction. Comparison of these strengths 
with the fully coupled collective strengths (solid curves) 
displays the collective El and E2 shifts. As a trivial 
fact, we observe the proper right shift for the isovec- 
tor GDR and the left shift for the isoscalar GQR. What 
is more interesting, the shifts for both resonances (in- 
cluding isovector DGR) depend on the isoscalar effective 
mass m*,/m. To explain this, one should remind that 
the smaller mj/m, the more stretched the single-particle 
spectrum (see also relevant examples in |5lj ) . And in- 



FIG. 2: The same as in Fig . 1 but for isoscalar GQR and 
(a, a') experimental data |50|. 



deed the El and E2 unperturbed strengths in the figures 
exhibit a systematic shift to higher energy while mov- 
ing from SkT6 to SkI3 (for El the shift is essentially 
weaker than for E2 Simultaneously, we have the 

corresponding evolution of the collective energy shifts. 
Namely, they decrease for GDR and increase for GQR 
(for exception of GDR for SkT6). All these trends lead 
to a remarkable net result: strong variations of the un- 
perturbed strengths and collective shifts with mj/ro con- 
siderably compensate each other so as the final GDR and 
GQR energies become much less sensitive to different 
Skyrme forces and approach the experimental values. 

The remaining energy differences of order ~ 1 MeV 
show some other, sometimes known, trends for the GR. 

First, we see a systematic (for exception of SkT6) 
downshift of the calculated GDR with increasing the 
symmetry energy a sym . This, at first glance, surprising 
result complies with experience of the systematic studies 
in spherical nuclei |5J] . The case of SkT6 looks as excep- 
tion from this simple rule. It has the same asymmetry 
energy as SkM* (see Table but a lower GDR reso- 
nance peak. This comes because the density dependence 
da sym /dp is abnormally low here (for reasons whose dis- 
cussion goes beyond the scope of the present paper). Be- 
sides the trend with a sym , Fig. n a l so hints a connection 
of the GDR energy with the sum rule enhancement fac- 
tor k and the related value of the isovector effective mass 
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m*/m. Namely, the smaller the effective mass (larger 
the k), the higher the GDR. The trends and connections 
mentioned above should, however, be considered with a 
bit of care. Indeed, the variation of a sym between the 
different Skyrme forces in Table [i] is rather small and 
probably not enough to demonstrate a strong and un- 
ambiguous trend. Besides, a complicated dependence of 
GDR on different isovector factors can spoil and entan- 
gle the concrete trends. In any case, analysis of the DGR 
trends deserves more systematic study which is now in 
reach with the efficient SRPA technique. 

Instead, the evolution of isoscalar GQR is more clear 
and systematic. We see a steady upshift of the GQR 
peak from SkT6* to SkI3, which complies with the known 
dependence on the isoscalar effective mass mj/m |55j| . 
namely, the lower the effective mass, the higher the 
GQR. This trend has a simple explanation. As was dis- 
cussed above, the low m^/m results in stretching the 
single-particle spectrum and thus in the upshift of the 
2qp quadrupole strength. In the GQR case, the op- 
posite dependence of the collective shift is not enough 
to compensate the strong 2qp upshift and hence we ob- 
tain the trend. Note that in our calculations SkT6 with 
rn^/m = 1 yields the best agreement for the GQR. This 
confirms to some extend the findings [55| for 208 Pb that 
a good reproduction of the GQR requires a large effective 
mass. 

Shape and width of GR in deformed nuclei are mainly 
determined by the deformation splitting and the spec- 
tral fragmentation due to interference with energetically 
close lph states (Landau fragmentation). Following Ta- 
ble UJ the different Skyrme forces provide quite similar 
quadrupole moments. Hence they result in close defor- 
mation splittings. Instead, the Landau fragmentation 
depends sensitively on the spectral pattern of a model, 
determined in a large extent by isoscalar and isovector 
effective masses. Our samples of forces contain sufficient 
variation of both ones. Nonetheless, it turns out that the 
width and shape of the GQR are practically the same for 
all the four Skyrme forces. 

At the same time, the strongest variation is found for 
the GDR whose width and gross-structure significantly 
depend on the force. It is seen that SkM* gives an arti- 
ficial third right peak (especially in 154 Sm) and thus an 
overestimation of the GDR width. This effect weakens for 
SLy6 (leading to the best description of GDR) and van- 
ishes at all for SkI3 (resulting in underestimation of the 
resonance width) . The appearance of the artificial right 
shoulder for SkM* and SLy6 was already noted for de- 
formed rare-earth and actinide nuclei Q and 208 Pb [Isj . 
For the particular Skyrme forces, this effect seems to be 
universal. It takes place for GDR in heavy nuclei, inde- 
pendently on their shape. As is shown in [37j, the right 
shoulder is provoked by an excessive collective shift and 
further enforced by the presence in the region of the 2qp 
bunch composed from the particular high-moment config- 
urations (7iT<79/2 and fl/iii/2 f° r 154 Sm and irlhn^ and 
iTi 13/2 for 238 U and 254 No) These configurations 



represent the intruder 1 + 1/2 states entering the valence 
shell due to the strong spin-orbital splitting. They form a 
dense 2qp bunch which is easily excited and thus provides 
high sensitivity of the right GDR flank. This feature can 
be useful for additional selection of the parameters of 
Skyrme forces. 

It would be very interesting to relate the above re- 
sults for the GDR width and profile with the effective 
masses. Some possible correlations can immediately be 
noted. For example, the evolution of the right shoulder 
from SkM* to SkI3 seems to correlate with a system- 
atic decrease of mj/m. Besides, the largest right shoul- 
der in SkM* can be related with the smallest value of 
m\/m for this force. This conclusion, however, requires 
still confirmation from more extended studies with an 
even broader basis of forces. At the same time, the re- 
sults hint that the properties of the isovector GDR prob- 
ably depend not only on the isovector m\/m but also on 
the isoscalar rn^jm. For the first glance, this statement 
looks surprising. However, we should take into account 
that ttiq/to influences the mean level spacings e in the 
mean field and hence can in principle affect the GDR. 
The resonance should depend on the ground state prop- 
erties which in turn are related to iri^/m. It worth also 
noting that both m^/m and m\/m are generated by one 
and the same term of the Skyrme functional (~ &i,6'i) 
and so probably are not fully decoupled in the dynamics. 
Some non-trivial relations between isoscalar and isovec- 
tor parameters are already discussed in literature, e.g. 
the relation Zra^jm ss Km\jm and its connection with 
cLsym |56| . And we see here an interesting field for future 
investigations. 

As a next step, let's consider contributions of the time- 
odd densities to the GDR and GQR. These results are 
illustrated in Fig. 3. The calculations show that only 
the current-current contribution (|52|) is essential while 
the contribution l|53|) connected with the spin density is 
negligible. So, in Fig. 3 we display only effect of the 
current density. First of all, it is worth noting that the 
most significant changes appear again at the right flank 
of the resonances. Probably, the high-moment I + 1/2- 
configurations mentioned above play important role not 
only for GDR but for GQR as well. 

The impact of the time-odd current is different in GDR 
and GQR. In the quadrupole resonance, we see the sys- 
tematic downshift and narrowing the strength. There is 
a clear correlation with the value of the isoscalar effective 
mass: the lower TOq/tti, the stronger the time-odd effect. 
After inclusion of the time-odd coupling, the GQR for 
different mj/m become much closer. So, the time-odd 
contribution weakens the influence of mj/m. 

In the dipole resonance, the time-odd shift is not so 
systematic. We observe no shift for SkT6, the upshift for 
SkM* and the downshifts for SLy6 and SkI3. Again one 
may note some correlation with effective masses. The 
specific SkM* case can be connected with very small 
m*/m for this force. Besides, one may note (for ex- 
ception of SkM*) increasing the time-odd impact with 
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FIG. 3: The isovector GDR and isoscalar GQR in 238 U cal- 
culated with the Skyrme forces SkT6, SkM*, SLy6 and SkI3. 
The strength is calculated with (solid curve) and without 
(dash curve) contribution of the time-odd current. 



lowering m^/m. In principle, the correlations between 
the influence of time-odd terms and effective masses in 
both GDR and GQR cases were more or less expected be- 
cause the current density enters the term of the Skyrme 
functional ~ b%, b[ just responsible for generation of the 
effective masses. This explains why the SkT6 case with 
mj/m = m*/m = 1 (no effective mass effects) does not 
demonstrate any time-odd impact. 

It is worth noting that the dominant contributions to 
the collective response from the principle terms of the 
Skyrme functional have different signs and thus, in a 
large extent, compensate each other (this can be easily 
checked in the SRPA by estimation of the different con- 
tributions to the inverse strength constants (|23[) . As a 
result, the smaller contributions (time-odd, spin-orbital 
and Coulomb) become important. 

The present study involves three nuclei from different 
mass regions and four Skyrme forces with various bulk 
properties. We have yet to disentangle more carefully 
the separate influences of them. This requires systematic 
variations of forces in a large set of test cases. Such 
systematic studies can be now easily performed due to 
the computationally efficient SRPA method. 



VI. CONCLUSIONS 

A general procedure for the self-consistent factoriza- 
tion of the residual nuclear interaction is proposed for 
arbitrary density- and current-dependent functionals. 
Following this procedure, the separable RPA (SRPA) 
method is derived. SRPA dramatically simplifies the cal- 
culations while providing a reliable description of nuclear 
excitations. The reduction of the computational effort is 
especially useful for deformed nuclei. In the present pa- 
per, SRPA with Skyrme forces is specified for description 
of collective dynamics in axially-deformed nuclei. 

For the first explorations, SRPA is applied for de- 
scription of isovector giant dipole resonances (GDR) 
and isoscalar giant quadrupole resonances (GQR) in de- 
formed nuclei from rare-earth ( 154 Sm), actinide ( 238 U), 
and superheavy ( 254 No) regions. Four Skyrme forces 
(SkT6, SkM*, SLy6, and SkI3) with essentially differ- 
ent bulk properties are used. The calculations show 
that SRPA can successfully describe multipole giant res- 
onances. A good agreement with available experimental 
data is achieved, especially with SLy6 for GDR. We did 
not find any peculiarities of the GR in superheavy nuclei. 
The behavior of the resonances in all three mass regions 
looks quite similar. 

We analyzed dependence of GDR and GQR descrip- 
tions on various Skyrme forces, in particular on the 
isoscalar and isovector effective masses and the symme- 
try energy. The contribution of the time-odd couplings 
was also explored. Some known trends for GDR with 
a sym and GQR with rrtg/m were reproduced. Besides, 
the close relation between the time-odd contribution and 
itiq jm was demonstrated for GQR. The calculations also 
hint some interesting (though not enough systematic) 
trends for GDR. Altogether, the results point out corre- 
lations between isoscalar and isovector masses and time- 
odd contributions. The correlation seem to be natural 
since all the items originate from one and the same term 
of the Skyrme functional. 

The time-odd and effective mass impacts manifest 
themselves mainly at the right flanks of the strength dis- 
tributions. The impacts are much stronger and diverse 
for GDR. Besides, this resonance exhibits the strong Lan- 
dau fragmentation. The GDR gross structure consider- 
ably depends on the applied Skyrme force and the related 
effective masses. As a result, the GDR structure can 
serve as an additional test for selection of the Skyrme 
parameters and related nuclear matter values. 
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APPENDIX A: DENSITY OPERATORS FOR 
SKYRME FUNCTIONAL 

In our study the Skyrme forces include time-even (spa- 
tial, kinetic-energy, spin-orbit) and time-odd (current, 
spin) densities associated with the hermitian operators 

^2 s (n - r), 

1=1 

V<5(r--r)V, 

»=i 

N a 

^2S(fl - r)Vx?, 
»=i 

1 N " 

^{v,S(fi-f)}, 



Ps(r 
T s {r 

Js(r 



where a is the Pauli matrix and N s is the number of 
protons or neutrons. 
The densities read as 



(Al) 



where Jf is the density operator, tpj{r) is the wave func- 
tion of the single-particle state j, and ig is the pairing 
occupation weight. 

APPENDIX B: WAVE FUNCTION IN 
CYLINDRICAL COORDINATES 

Cylindrical coordinates p, z, 8 are defined as 

x = pcosd, y = psind, z = z . 

Then the single-particle particle wave function and its 
time reversal have the form of spinors 



(Bl) 



tpj(r) = TipAr) = [ H J ^> z > e ( B 2) 

where Kj is the projection of the complete single-particle 
moment onto symmetry z-axis of the axial nucleus. 

Expressions for differential operators in cylindrical co- 
ordinates are elsewhere, see e.g. |46| . 

APPENDIX C: PAIRING CONTRIBUTION 



The pairing is treated with Bardin-Coopper-Schriffcr 
(BCS) method. Then the SRPA values gain the pairing 



factors involving coefficients Vj and Uj of the Bogoliubov 
transformation from particles to quasiparticles. 

The densities and currents in the Skyrme functional 
dominate in the particle-hole channel. The pairing den- 
sity falls into another channel provided by the two- 
particle excitations. It involves a different pairing weight, 
namely UjVj rather then the vq for the standard densities 
I|A1(I . Specifically the pairing density reads 



(CI) 



In the case of pairing, the hermitian one-body opera- 
tors with a given time-parity ((7™ = 1,7t = — 1) obtain 
in the p/i-channel the form 

A = 2 y £<ij\A)u# ) (At j +A ij ), (C2) 

ij 

B = 2^2(ij\B\)u£\At+A ij ), (C3) 



where 



Ah — a] at, AiA = cbcb 



are two-quasiparticle operators and 



V,\j = UiVj + UjVi, V,\j ) = UiVj - UjVi 



(C4) 



(C5) 



are the pairing factors. This is the case for time-even 
operators Q s k and X s k and the time-odd operator Y s k- 
The time-odd operator 

P sk = i[H, Q sk ] = i{[ho, Q sk ] + [K7, Qsk]} 

= i[h ,Qsk]~Y s k (C6) 

have a more complicated form because of the additional 
term i\ho,Q s k\- Namely it reads 

Psk = 2^{i2e ijU { +\ij\Q sk )-utf{ij]Y: k )} 

ijes 



(C7) 



The SRPA formalism in sections ITTland lTTTl is presented 
in a general form equally suitable for the cases with and 
without pairing. In the case of pairing, the ground and 
perturbed lph many-body wave functions are replaced 
by their BCS counterparts, the summation indices p and 
h run the quasiparticle states, and the involved values 
(densities, operators and matrix elements) acquire the 
pairing factors given above. 

The pairing is not important for giant resonances con- 
sidered in the present study. So, we freeze the pairing in 
the dynamics and do not present here the explicit form of 
the pairing contribution (pairing vibrations) to the resid- 
ual interaction. Some examples of this contribution can 
be found in |22|. 
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